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Abstract 



The effect of stellar scintillation on the accuracy of photometric measurements is ana- 

r^ • lyzed. We obtain a convenient form of estimaton of this effect in the long exposure regime, 

when the turbulence shift produced by the wind is much larger than the apert ure of the 



Q ■ telescope. A simple method is proposed to determine index S3 introduced by Ke nvon et alJ 



( 20061 ). directly from the measurements with the Multi Aperture Scintillation Sensor (MASS) 



^ . without information on vertical profile of the wind. The statistics S3 resulting from our cam- 

paign of 2005 - 2007 at Maidanak observatory is presented. It is shown that these data can 
be used to estimate high-altitude winds at pressure level 70 - 100 mbar. Comparison with 
the wind speed retrieved from the NCEP/NCAR global models shows a good agreement. 
Some prospects for retrieval of the wind speed profile from the MASS measurements are 
outlined. 
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1 Introduction 

The accuracy of measurements of brightness of astronomical objects is affected by many factors. 
It is clear that they are influenced in different ways depending on the applied photometric 
5^ \ techniques. One of these effects is the stellar scintillation caused by optical turbulence (OT) 

in the Earth's atmosphere. Unlike other variations of the atmospheric transparency, the light 
fluctuations due to the scintillation do not correlate for the images separated by more than the 
isoplanatic angle (2" — 5" in the optical range). Therefore, errors introduced by this effect appear 
always and increase when the brightness of two stars is c ompared. 

This problem was noted and studied a long time ago ( Younel . Il969l : iDravins et al.l . Il997a| ]bl. 



19981 ). Recently, interest in this problem manifested itself again, but in the overall context of a 



study and ch aracterization of the OT in various astronomical observatories and prospective sites. 
In the paper iKenyon et al.l ( 20061 ) a special index ^3 have been introduced to relate the vertical 



distribution of the turbulence and wind speed to the variations of the light flux as measured in 
a large telescope. 

Monitoring of this parameter and its analysis is very useful for comparing different sites and 
for predicting or planning of photometric observations. However, the approach used in cited 
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work i mplies prior restoration of the OT profi l es fro m the MASS - Multi Aperture Scintillation 
Sensor ( Kornilov et al.1 . |2003J; iTokovinin et alj . |2003| ) data and the use of additional information 
on wind speed. This is a very intensive procedure, besides burdened by additional errors 

In this paper we propose a method to determine the index ^3 directly from the MASS 
measurements without using additional information. The presented method is preceded by a 
theoretical consideration of the scintillation in the temporal domain. 



2 Temporal averaging of the scintillation 

The approximation of weak perturbations implies the independence of the wavefront distortions 
produced in a turbulent layer from the turbulence located above it. This leads to the fact that 
the scintillation index s^ — the variance of relative fluctuations of the light intensity, is described 
by the sum of the scintillation indices produced by individual layers: 



Clih)Wih)dh 



il) 



where W{h) is the weighting function (WF) which depends on the size and shape of the receiving 
aperture and does not depend on the altitude distribution of the structural coefficient of the 
refractive index C^(/i). WF represents a power of the scintillation generated by a layer of unit 
intensity located at a height h. 

The MASS method involves simultaneous measurement of the scintillation indices in 4 con- 
centric apertures of different diameters, leading to 4 normal and 6 diffential scintillation indices. 
The v ertical OT profile is restored from the mea sured indices and the theoretically calculated 
WFs ( Tokovinin et al.l . l2003l : lKornilov et al.l . 120071 ). The calculation of the set of functions W{h) 
assumes that the light intensity measurement has a "zero exposure", i.e. the averaging factor 
is related only to the receiving aperture. This is not quite true in real measureinents, so the 
problem of finite time of the measurement was considered separately by iTokovininI (120021 ) . 

The question of t he quantitative effect of averaging on the scintillation was considered in 



many papers, see e.g. iDravins et al.l (Il997al). We consider below this pro blem using the termi 



no logy of the MA SS theorv (JTokovinin et al.l . l2003l : iKornilov et al.l . 120071 ) and the results given 
bv lMartinI ( 19871 ) and iTokovinirj ( 20021 ) in order to obtain formulas more convenient for practical 
application. 

It is obvious that the variance o"^ of the relative fluctuations of the light intensity measured 
with finite exposure time is determined by the integrated effect of all turbulent layers on the 
line of sight, as in the case of zero exposure ([1]), but the WFs will be different. 

Assuming t hat th e temporal evolution of the distorted wavefront is described by its transla- 
tion (Taylor's ( 19381 ) hypothesis of the frozen turbulence) by wt, we can describe the effect of 
temporal avera ging by a convolution of the re ceiving aperture with a line of wind shift in the 
spatial domain ( Tokovinirj . 120021 : iMartirj . (19871) . In calculating the WF as integral of the filtered 
spatial power spectrum of scintillation, an additional multiplication by the temporal averaging 
filter appears. 

The new WFs W'{w,T,h) will depend not only on the layer altitude h, as previously, but 
also on the wind speed w = w{h) and averaging time r. 



W'{w,T,h) = 9.62 A" 



d/ /-8/3 sin2(7rA/i/2) A{f) As{w, T, /), 



where / is the modulus of spatial frequency, A{f) is the aperture filter (axi-symmetric for the 
MASS apertures) and As{w,t, f) is the additional spectral filter of the wind shear along the 



axis X, which is the square of the Fourier Transform of the rectangular window of the size wt. 
In polar coordinates it can be written as 



Asiw,Tj) 
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)sinc {frw cos 4>) =Ti{fTw), 



(3) 



where the function 71 (^) can be expressed through the Bessel functions Jq and_Jjand_the 
Struve funct i ons H n and Hi (the definition of Struve functions can be found in (jStruvd . Il882l : 
Abramowita . 1 19651 )): 



ri(C) = 2Jo(27rO - ^^^^ - 7rJo(27rOi^i(27re) + 7rJi(27rOi^o(2vre) 



(4) 



The behavior of this function is shown in Fig. [TJ Although the paper of lTokovininI (J2002l ) does 
not provide the analytical expression of this function, its correct asymptotic was given. We only 
refine the interval of the asymptotic. For small S,, the function 71 (^) ~ I — tt'^S.'^/Q, as follows from 
its series expansion in the neighborhood of 0. The quadratic approximation provides accuracy 
better than 0.02 until the vr,^ < 1. When ^ — )• oo, the function 71(^) ~ I/tt^, and starting from 
^ ~ 1 the relative difference is less than 0.04. 
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Figure 1: Wind shear spectral filter As{w,T,f) = Ti{fTw) 

The asymptotic behavior of the function at ^ — )• and ^ — )• oo can be used to describe the 
two limiting cases: the short exposure an d long exposure r egimes. The regime of short exposure 
Tw <C D have been analyzed in detail by [Tokovininj ( 20021 ). 



3 Long exposure regime 

Classical astronomical photometry almost always works in the approximation of the long expo- 
sure TW ^> D. In this situation, frw ^ 1 and one can use the asymptotics 71 (^) ~ I/tt^. Then 
the WF looks like 



W'{w,T,h) =9.62 A" 
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where U'{h) is the value of the integral over spatial frequency. Note that the integrand has 
no singularities and is integrable without any problems, as well as in the zero-exposure case. 
If / —7- 0, it tends to 0, justifying the use of this approximation. It is evident that the WF 
W'{w, T, h) has a simple structure and the scintillation power averaged over exposure r is defined 
by the following expression: 

7o w{h)T 

In terms of MASS the method, this expression has a very convenient form that allows, in 
principle, to determine the vertical distribution of C'^{h)/w{h) from measurements made in a 
long-exposure regime. Note that the approximation of long exposure is valid at lower values of 
r for small apertures, such as in the MASS device. 

In the limiting case of infinitely small aperture U'j h) = 13.52 A~^/ ^ h'^'^ can be obtained 
using the fact that A{f) = 1 and changing the variables (JRoddieri . ll98ll ). For the opposite case 



of large aperture D ^ rp (Fresnel scale rp = vA/i), the function [7^ (/t) can be calculated by 
replacing the Fresnel filter sin^(7rA/i/^) with (ttA/i/^)^, as done by iRoddien ( 198ll ). Turning to 



the dimensionless frequency q = fD, we find that U'{h) = 10.66 D ^'^/i^ or 

a\T) = WmD-''/\-' r^M^^h (7) 

Jq w{h) 



If -D = 1 and t = 1, this expression represents the index S? introduced in paper iKenvon et al. 



(|2nn6l ). 



In the case of long exposures, the dependence of scintillation on airmass M^ (or zenith 
angle) is obvious: o"'^(r) oc M^. However, in the formulas ([6]) and ([7]) w is the component of the 
wind velocity perpendicular to the line of sight. Naturally, the wind speed is normally directed 
horizontally. If the azimuth of the wind differs from the azimuth of the line of sight by ±90°, 
then the perpendicular component is simply equal to the absolute speed. If the wind direction 
coincides with the plane of the line of sight, then perpendicular component is w_\_ = w cos z and 
the dependence of o"^(r) on the air mass increase s to M^. The dependence of the scintillation 
power on azimuth was described bv lYouneJ ( 19691 ). 



4 Direct way to calculate the 6*3 index 

The value 51 is the variance of the stellar scintillation measured at 1 m telescope with 1 s 
exposure. The main point of the method for calculating the index 5*3 (needed for estimating the 
photometric accuracy on large telescopes) directly from the MASS data is that S"! is represented 
as a linear combination of long-exposure scintillation indices at different apertures of the MASS 
device. Unlike the instantaneous indices s^, they will be denoted as <t^. 

In the MASS data processing, the so-called atmospheric moments are calculated. In partic- 
ular, the second moment is 

poo 

M2 = / Cl{h) h^ d/i, (8) 

The problem of approximation of the h? dependence is solved in the MASS data processing 
by a linear combination of the WFs W{h). The case under study here differs only in that the 
approximation of h? must be done with a set of U'{h) functions, having the dimension of m?'^ . 
A set of U'{h) is shown in Fig. [2l These functions are similar to the usual WFs, but are more 
smooth owing to the suppression of high spatial frequencies. It is clear that the time-averaged 
scintillation are much smaller. Therefore, the values of the functions U' and W differ by about 
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Figure 2: MASS WFs W{K) (dashed) and \J'(K) (solid lines). Steady curves represent normal 
indices, saturated ones — differential indices. Dot-dashed line shows a square-law altitude 
dependence 

2 orders of magnitude. Both sets of WFs are computed for the polychromatic case with the 
spectral curve of the original MASS device and an AO V star. 

The curve for the A aperture differs slightly from the WF of an infinitely small aperture at 
A = 0.5/im (effective wavelength). Strictly speaking, the limiting behavior for large apertures 
is also valid for the monochromatic case. For example, the real U'{h) for the photometric band 
V grows with h more slowly. Moreover, telescope central obscuration also slows the growth of 
U'{h). However, unlike the case of short exposure, these effects are small. Even for a 0.3-m 
aperture, the ratio of U'{h) to 10.66 D~^'^h'^ differs from 1 by no more than 5% at 12 km altitude 
for the R band and with the central obscuration of 33%. In the case of a 0.5-m telescope, all 
the differences do not exceed 2% at all altitudes. 



5 Approximation of the quadratic dependence with the MASS 
WFs. 

For a rough estimate of the S3 index one can use an evaluation of o"^ from measurements in 
the D aperture of the MASS device. Although its diameter is about 11 cm and the central 
obscuration is more than 50%, the corresponding WF aproximates the quadratic dependence 
quite well. The maximum difference is achieved at high altitudes (see Fig. [3]) and is f«25% at 
the tropopause altitude. This means that the ^3 error will not exceed ~15%. 

The accuracy of the quadratic approximation can be significantly improved by representing 
the WF of a 1 m telescope as a linear combination of functions U'{h). In Tabled] the coefficients 
of this expansion are listed. The coefficients Cj are obtained by the decomposition on the whole 
set of U'{h) functions, dj — only on the 4 functions corresponding to normal indices in apertures 
A, B, C and D. Both approximations are shown in Fig. [3l where the ratios of the approximated 
function to 10.66 D~'^'^h'^ are plotted. It is clear that both approximations provide a similar 
accuracy, better than 5%, over the entire range of altitudes 0.5 - 30 km, and better than 2% in 



Table 1: Approximation of h'^ by the WFs U'{h) of the original MASS instrument. 
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Figure 3: The ratio of the approximating linear combinations of the WFs to 10.66 /i^. Black 
dots — the approximation using only normal indices, empty circles — using all 10 indices. The 
dashed curve is the ratio of the WF for the D aperture of the device (11 cm) 



the tropopause region. For the MASS/DIMM instrument the situation is slightly worse because 
the typical size of the aperture D is 8 cm and the central obscuration is larger. 

In view of the linear relation between the indices and their corresponding WFs d?]) the 
following formula is correct: 
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(9) 



and the calculation the S^ index reduces to the measurement of the variance of the flux in the 
MASS apertures averaged over 1 s. 



6 Calculation of the S^ index 

The required quantities are not calculated by the MASS software during measurements. How- 
ever, all the necessary information is stored in the output files * . stm, where the intensities in 
apertures A, B, C and D averaged over 1-s exposure are recorded. Note that these values in 
the file are given for the 1-ms exposure, in fact the counts are 10^ -i- 10^ and, consequently, the 
contribution of the photon noise is small. 
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Figure 4: Diffential distributions of the varinces of 1 s average fluxes in the A, B, C and D 
apertures of the MASS instrument 

Figure 5: Distribution of the 5*3 index. The sohd hue is obtained for variances estimated over 1 
minute, dashed hue — over 4 minutes. The stair-Uke Hne is the diffential distribution of S3 



Unfortunately, for a typical 1 minute acquisition time there are only 60 measurements. 
Therefore, the relative accuracy of the computed variance is approximately 25%. For more 
accurate values, further averaging (over 4-minute interval) is needed. 

Selection of valid data presents some problems because the fluctuations on 1-s temporal scale 
can include significant intensity changes of another origin — due to changes in transparency, 
clouds, or large guiding errors. 

To investigate the atmospheric limit of photometric accuracy, the data of the measurements 
performed on the mo untain Maidanak in 2005 - 2007 with the MASS instrument were used 
( Kornilov et al.l . l2009l ). The variances a'^ for all apertures were calculated with help of the mean 
square difference of adjacent 1-s means 
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where F is the mean of the 1-s intensities {Fi} in a given aperture. The last term is the 
contribution of the photon noise. 

Use of the differences effectively suppresses the contribution of intensity variations on time 
scales longer than 1 sec, of non-scintillation nature. At these scale, the scintillation itself is 
uncorrelated, which leads to the asymptotics ct^ oc t~^ . We used only data where the average 
count in the channel D was greater than 100 pulses/ms and the relative variance did not exceed 
0.002. 

Differential distribution of thus defined scintillation indices for the entire campaign period 
(45 714 1-min points, or about 80% of all measurements) are shown in Fig. SI Values are reduced 
to the zenith. Medians and quartiles of the distributions of 1-s indices in the apertures A, B, C 
and D are given in Table. [2j The shape of the differential distributions is as expected and shows 
that data censoring does not distort them. 



Table 2: Characteristics of the distributions of the variances 1 s mean fluxes and the indices 5| 
and 53 



A 



B 



C 



D 
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"^3 



^3 



quartile 25% 


0.000246 


0.000208 


0.000152 


0.000088 


5.83-10-^ 


0.00242 


median 


0.000386 


0.000322 


0.000234 


0.000135 


9.07-10-*^ 


0.00302 


quartile 75% 


0.000589 


0.000488 


0.000351 


0.000201 


13.65-10"*^ 


0.00370 



The 5*1 index was calculated using the coefficients dj from Table. [TJ The median value 
)-5 ];^4/3g_ fjij^g index computed over 4- minute intervals has a slightly higher 



of Si is 0.91 



median of 0.93-10^ m ''^s. The difference is negligible, proving that the accuracy of the variance 
estimation over 1 minute accumulation time is sufficient. Clearly, this accuracy is mainly defined 
by the error of the variance in the D aperture. Cumulative distribution of the index ^3 is shown 
in Fig. [5j Naturally, this distribution is narrower than the distribution of S^. The quartiles of 
this distribution are listed in Table. [2l 

The median of the ^3 distribution equals 0.0030 m^'^s^'^. This value is in good agreement 
with the estimation of 53 for the Cerro Pachon and Cerro Tololo observatories from the paper 
Kenvon et al.l ( 20061 ). where the medians of 0.0028 and 0.0032 m^/'^s^'^, respectively, are given. 

The evolution of the quantity 5*3 during the 2005 - 2007 campaign is shown in Fig. [U 
where the medians over each ~250 min data points are plotted. The annual cycle is clearly 
visible. Seasonal variations are considerable, the index 5*3 varies by almost two times. The best 
photometric accuracy is expected in the period from October to April. 
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Figure 6: Seasonal behavior of the median values of the index 5*3 



7 Evaluation of the effective wind speed 

If we compare the formulas ^ and ([8|) , we can see that their ratio gives some kind of estimation 
of the wind speed w. More precisely, the ratio of the second atmospheric moment M2 to the 
index S^ gives the average value of l/w{h) weighted by C^{h) h?: 

1 \"^ 10.66 -Mo 



The M2 is calculated in the sta ndard MASS data p rocessing and, apart from a coefficient, 
equals the index S^ introduced in ( Kenvon et al.l . I2OO6I ) . The weight C^{h) h? reaches maximum 



usually at altitudes above 16 km and the expression (fTTIl evaluates thus high-altitude wind. 
Additional uncertainty is associated with averaging of reciprocals, when the resulting average 
may be underestimated. At Maidanak, the height of 16 km above the summit corresponds to 
the barometric altitude of 18 - 19 km, or the pressure level of ~ 70 mbar. 

The result of high-altitude wind speed calculation are shown in Fig. [71 The distributions 
of the wind speed vjq and f 100 at pressure levels of 70 and 100 mbar as retrieved from the 
database NCEP/NCAEo for moments of the measurements are also plotted. It is seen that 
the distribution of w is bracketed by these two distributions, indicating the effective height 
determined by the C^{h) h? weighting is indeed somwhere between these two levels. 

Comparison of the individual points of w calculated for the values of 5"! and M2 estimated 
over 4-min. with vjq interpolated at these moments is shown in Fig. [8] in the form of a plot 
of the median values over groups of 250 points. A significant correlation and an almost linear 
relation between w and W70 are clearly visible. In the weak wind domain a systematic excess of 
w over ^70 is observed. This may be due to the fact that in such situations the contribution of 
lower layers with faster wind and stronger turbulence becomes important. 

In Fig. [9] the long-term behavior the of high-altitude wind speed above the mountain Maid- 
anak is shown. It is seen that the minimum value of high-altitude wind w ~ 10 m/s is observed 
from early July to mid-October. A similar behavior is characteristic of W70, but not of wind at 
the pressure level of 200 mbar. 

A combination of the WFs U'{h) that grows with h slower than h? (e.g. linearly) can be 
found. However, since the WFs of the normal scintillation indices increase with altitude faster 
than /i^'^, such approximation causes a noticeable increase in the noise. It is preferable to use 
the temporal index in the A aperture (oc /i^'^) and the atmospheric moment M4/3, computed 
from the indices s^. In any case, the dominant layer will be located lower than in the case 
considered above and the ratio similar to the formula (llip is likely to reflect the wind speed in 
the tropopause. 

8 Conclusion 

A method to obtain statistical prediction of photometric accuracy based only on the data from 
the MASS instrument is proposed. Cor uparison of the resul ts obtained by this method and by 
integration of the OT and wind profiles ( Kenyon et al.l . 120061 ) can be used to verify wind models 



above astronomical sites. 

It is particularly interesting to perform such a comparison for the measurenients w ith MASS 
instrument performed in Antarctica at Dome C, as described in ( Kenvon et al.l . 120061 ). There 



IS 



some concern that their data censoring procedure, leaving only 16% of useful data, may bias the 
53 estimates. An ability to estimate wind speed at pressure level of 70 - 100 mbar directly from 



^http://www.esrl. noaa.gov/psd/data/reanalysis/reanalysis.shtml 
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Figure 7: Cumulative distributions of the winds. Solid line — wind evaluated w by formula 
PT|) . dash-dotted — wind at 100 mbar pressure level, dashed — at 70 mbar level selected for 
moments of the measurements. Stair line — differential distribution of the w 

Figure 8: Relation between the w and wind at 70 mbar pressure level. Empty circles are meadian 
values in 250 points groups 
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Figure 9: Seasonal behavior of the wind w (black dots) and wind at 70 mbar W70 (empty circles) 
for the entire period of measurements with the MASS instrument on Maidanak 
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the MASS data will help to check the applicability of the wind models and will avoid biases 
caused by replacing instantaneous values by their statistical estimates. 

Additionally, we note here that 1-s averaging time in the case of MASS apertures is too long. 
A regime of the long exposure with the relevant asymptotics starts at shorter exposures on the 
order of 0.1 — 0.2 s, and the use of such short time would markedly improve the accuracy and 
reliability of the temporal index. The behavior of WFs U'{h) (Fig. [2]) indicates that by using 
a set of precise measurements of temporal indices and the restored C^{h) profile, a reasonable 
solution of the system equations like ^ for the unknowns l/w{h) may be achieved. 

Realization of this possibility will permit 1) to obtain additional information about the 
current state of the turbulent atmosphere, and 2) to correct more precisely the MASS/DIMM 
data for the effects of finite exposure time. 

The author thanks his colleagues who participated in the campaign to measure the optical 
turbulence above the mountain Maidanak in 2005 - 2007, producing the data used here. The 
work was partially supported by the RFBR grant 06-02- 16902a. 
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